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Abstract 

We  study  the  one  machine  scheduling  problem  with  release  and  delivery  times  and  the  mini¬ 
mum  makespan  objective,  in  the  presence  of  constraints  that  for  certain  pairs  of  jobs  require 
a  delay  between  the  completion  of  the  first  job  and  the  start  of  the  second  (delayed  prece¬ 
dence  constraints).  This  problem  arises  naturally  in  the  context  of  the  Shifting  Bottleneck 
Procedure  for  the  general  job  shop  scheduling  problem,  as  a  relaxation  of  the  latter,  tighter 
than  the  standard  one  machine  relaxation.  The  paper  first  highlights  the  difference  between 
the  two  relaxations  through  some  relevant  complexity  results.  Then  it  introduces  a  mod¬ 
ified  Longest  Tail  Heuristic  whose  analysis  identifies  those  situations  that  permit  efficient 
branching.  As  a  result,  an  optimization  algorithm  is  developed  whose  performance  is  com¬ 
parable  to  that  of  the  best  algorithms  for  the  standard  one  machine  problem.  Embedding 
this  algorithm  into  a  modified  version  of  the  Shifting  Bottleneck  Procedure  that  uses  the 
tighter  one  machine  relaxation  discussed  here  results  in  a  considerable  overall  improvement 
in  performance  on  all  classes  of  job  shop  scheduling  problems. 


1  Introduction 


The  problem  that  we  address  in  this  paper  arises  in  the  context  of  general  job  shop  scheduling. 
In  the  job  shop  scheduling  problem,  jobs  are  to  be  processed  on  machines  with  the  objective  of 
minimizing  some  function  of  the  completion  times  of  the  jobs,  subject  to  the  constraints  that 
(i)  the  sequence  of  machines  for  each  job  is  prescribed,  and  (ii)  each  machine  can  sequence 
at  most  one  job  at  a  time.  The  processing  of  a  job  on  a  machine  is  called  an  operation,  and 
its  duration  is  a  given  constant.  The  objective  chosen  here  is  that  of  minimizing  the  time 
needed  to  complete  all  the  jobs,  called  the  makespan. 

Let  N  =  {0,1,..., n}  denote  the  set  of  operations,  with  0  and  n  the  dummy  ‘‘start’" 
and  “finish”  operations,  respectively,  M  the  set  of  machines,  A  the  set  of  pairs  of  operations 
constrained  by  precedence  relations  representing  condition  (i)  above,  and  Ek  the  set  of  pairs 
of  operations  to  be  performed  on  machine  k  and  which  therefore  cannot  ov^erlap  in  time,  as 
specified  in  (ii).  Further,  let  dj  denote  the  (fixed)  duration  or  processing  time,  and  tj  the 
(variable)  start  time  of  operation  j.  The  problem  can  then  be  stated  as 

minimize 


subject  to 


tj-ti  >  d,  [ij)  e  A 

{P)  u  >  0  I  e  N 

tj  —  U  >  di  V  ti  —  tj  >  dj  {i,j)  G  Ek,  k  E  M. 

Any  feasible  solution  to  (P)  is  called  a  schedule. 

It  is  useful  to  represent  this  problem  on  a  disjunctive  graph  G  :=  {N,A.E)  with  node 
set  N,  (directed)  arc  set  A,  and  (undirected,  but  orientable)  edge  set  E.  The  length  of  an 
arc  {i,j)  G  A  is  di,  whereas  the  length  of  an  edge  {i,j)  G  P  is  either  </,  or  dj,  depending 
on  its  orientation.  Each  machine  k  corresponds  to  a  set  Nk  of  nodes  (operations)  and  a 

set  Ek  of  edges  that  together  form  a  disjunctive  clique.  Figure  1  shows  a  disjunctive  graph 

representing  a  job  shop  scheduling  problem  with  3  jobs  and  3  machines.  The  disjunctive 
cliques  corresponding  to  the  3  machines  have  node  sets  Ni  :=  {1,5.8},  !\j  :=  {2,4,9}, 
Ns  :=  {3,6,7},  with  their  edges  drawn  in  dotted  lines. 
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Figure  1;  Disjunctive  graph  for  a  3-job  3-machine  problem 

The  job  shop  scheduling  problem,  known  to  be  strongly  NP-hard.  is  among  the  hardest. 
This  is  dramatically  illustrated  by  the  fact  that  an  instance  involving  10  jobs  and  10  machines, 
proposed  by  Fisher  and  Thompson  [FT63]  in  1963,  remained  unsolved  for  more  than  a  quarter 
of  a  century,  although  every  available  algorithm  was  tried  on  it.  This  makes  it  imp-  rative 
to  design  efficient  approximation  methods. 

The  Shifting  Bottleneck  (SB)  Procedure  [ABZ88]  is  one  such  method.  It  sequences  the 
machines  consecutively,  one  at  a  time,  with  the  remaining  unsequenced  machines  ignored 
(i.e.  the  corresponding  edge  sets  removed)  and  the  machines  already  .sequenced  held  fixed 
(i.e.  the  corresponding  edges  replaced  by  directed  arcs).  At  each  step  a  bottleneck  machine  is 
determined  from  among  those  not  yet  sequenced,  by  sc>lving  a  one  machine  scheduling  prob¬ 
lem  for  each  unsequenced  machine  and  choosing  the  one  with  the  maximum  makespan.  This 
bottleneck  machine  is  then  sequenced  optimally  by  using  the  solution  to  the  corresponding 
one  machine  problem.  Once  all  the  machines  have  been  sequenced,  each  machine  in  turn  is 
freed  up  and  resequenced,  with  the  .sequences  on  the  remaining  machines  held  fixed. 

While  the  Shifting  Bottleiieck  Procedure,  having  been  extensively  tested,  was  found 
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more  efficient  than  any  version  of  the  traditional  dispatching  rule-based  heuristics,  whether 
deterministic  or  randomized,  several  researchers  (see  Dauziere- Peres  and  Lasserre  [DL90].  as 
well  as  Tiozzo  [TsS])  found  the  following  weakness  or  shortcoming  in  the  procedure.  When 
the  procedure  fixes  the  sequence  on  a  machine,  it  may  thereby’  create  a  precedence  constraint 
between  some  pair  of  jobs  on  some  unsequenced  machine.  In  terms  of  the  disjunctive  graph, 
orienting  the  edges  of  a  disjunctiv'e  clique  (i.e.  replacing  them  with  directed  arcs)  creates 
new  paths  in  the  graph,  some  of  which  may  join  two  nodes  of  a  disjunctive  clique.  If  a 
(directed)  path  is  created  between  nodes  i  and  j  of  a  di,sjunctive  cliciue,  this  implies  the 
constraint  tj  —  ti  >  L{i^j),  where  L{i,j)  is  the  length  of  the  path  from  i  to  j.  For  instance, 
if  in  the  graph  of  Figure  1  the  edges  joining  the  node  set  Ni  ;=  {1,5,8}  are  replaced,  as  a 
result  of  sequencing  machine  1,  by  the  directed  arcs  {(5, 1),  (5,  8),  (1, 8)},  that  replacement 
creates  the  path  {(4, 5),  (5, 1),  ( 1, 2)}  from  4  to  2,  of  length  L(4, 2)  =  +  d\.  which  in 

turn  imposes  the  constraint  1-2  ~  L{A,2)  (see  Figure  2). 


Figure  2:  Disjunctive  graph  after  sequencing  machine  1 

Note  that  the  constraint  L{i,j)  is  different  from  the  standard  precedence  con¬ 

straint  tj  —  ti  >  d,  that  would  result  from  arbitrating  the  disjunction  fj  — >  d,  V  t,  —  fj  >  d,. 
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in  that  the  delay  L{i,j)  involved  in  the  constraint  is  different  from  (typically  larger  than) 
di,  the  delay  involved  in  the  standard  precedence  constraint.  To  make  this  distinction,  we 
call  the  constraints  described  here  precedence  constraints  with  delay  or,  for  short,  delayed 
precedence  constraints  (DPC). 

As  [T88]  and  [DL90]  correctly  observed,  the  appearance  of  DPC’s  during  the  application 
of  the  SB  Procedure  vitiates  the  intent  of  the  heuristic  and  may  negatively  affect  the  outcome. 
Indeed,  by  ignoring  DPC’s  while  solving  a  one  machine  problem,  we  are  solving  a  less 
constrained  problem  than  we  should,  hence  we  may  make  an  unintended  choice  for  the 
bottleneck  machine  and  for  the  job  sequence  on  that  machine;  all  of  which  may  negatively 
affect  the  quality  of  the  schedule  found  by  the  SB  Procedure.  To  remedy  this  situation, 
Dauziere- Peres  and  Lasserre  [DL90]  developed  a  heuristic  for  the  one  machine  problem  with 
DPC’s.  Using  this  heuristic  instead  of  Carlier’s  exact  algorithm  for  solving  the  one  machine 
problems  that  arise  in  the  SB  Procedure,  in  other  words  substituting  an  approximation 
method  applied  to  a  better  chosen  problem  for  the  exact  method  applied  to  the  earlier 
problem,  they  managed  to  obtain  better  overall  results. 

Our  purpose  in  this  paper  is  to  develop  an  optimization  algorithm  for  solving  the  one 
machine  problem  with  DPC’s,  in  order  to  further  improve  the  SB  Procedure.  Of  course, 
using  an  optimization  algorithm  only  pays  if  the  gain  in  the  quality  of  solutions  to  the 
one  machine  problem  with  DPC’s  is  not  outweighed  by  the  excess  of  computational  effort 
required.  The  kind  of  one  machine  problem  solved  in  the  SB  Procedure  (i.e.  with  release 
and  delivery  times)  is  strongly  NP-complete  even  without  DPC’s;  nevertheless,  it  can  be 
solved  in  most  cases  quite  efficiently  due  to  the  possibility  of  using  a  branching  rule  that 
skips  large  subsets  of  potential  search  tree  nodes  (see  Potts  [P80]  and  Carlier  [C82]).  That 
same  efficient  branching  rule  can  also  be  used  if  some  standard  precedence  constraints  (as 
opposed  to  DPC’s)  are  imposed,  i.e.  if  some  of  the  disjunctive  arc  pairs  are  arbitrated  before 
.solving  the  problem.  However,  this  is  no  longer  the  case  in  the  presence  of  DPC's. 

The  main  accomplishment  of  this  paper  is  an  analysis  of  the  structure  of  the  one  machine 
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scheduling  problem  with  DPC’s  that  identifies  the  situations  in  which  the  above  mentioned 
“strong”  branching  rule  can  legitimately  be  used.  As  a  result,  an  algorithm  is  developed 
which  uses  the  strong  branching  rule  whenever  possible,  and  an  alternativ'e  rule  when  nec¬ 
essary;  with  the  practical  outcome  that  the  overall  efficiency  of  the  algorithm  is  comparable 
to  that  for  the  standard  one  machine  problem. 

The  rest  of  the  paper  is  organized  as  follows.  The  next  section  (2)  discusses  the  relation 
between  the  standard  one  machine  problem  and  its  counterpart  involving  delayed  precedence 
constraints,  highlighting  the  difference  through  some  relevant  complexity  results.  Section  3 
states  and  analyzes  the  Longest  Tail  Heuristic  for  the  one  machine  problem  with  DPC’s. 
Its  focus  is  on  identifying  the  situations  in  which  the  Pott.s-Carlier  property  still  holds,  thus 
permitting  efficient  branching.  Section  4  describes  a  branch  and  bound  algorithm  based  on 
the  foregoing  analysis.  Finally,  the  last  section  (5)  discusses  computational  results  with  the 
algorithm  both  as  a  stand-alone  and  as  part  of  an  improved  SB  Procedure. 

Some  of  the  results  of  this  paper  were  presented  in  [BV91]. 

2  The  One  Machine  Problem  with  DPC’s  Versus  the 
Standard  One  Machine  Problem 


At  some  point  during  the  SB  Procedure,  let  Mo  be  the  set  of  machines  already  sequenced, 
and  let  D{Mo)  be  the  directed  graph  obtained  from  G  by  replacing  all  disjunctive  arc  sets 
corresponding  to  the  machines  in  Mo  with  the  conjunctive  arc  sets  representing  the  sequences 
chosen  for  those  machines,  and  deleting  all  disjunctive  arc  .sets  corresponding  to  the  machines 
in  M  \  (Mo  U  {A;})  for  some  fixed  k  ^  M  \  Mo.  Then  the  standard  one  machine  scheduling 
problem  to  be  solved  for  machine  k  can  be  stated  as 

minimize 

subject  to 

tn  —  ti  ^  d,  +  qi  i  I 

ti  >  r,  i  €  / 

ij-t,>d,  V  t,-tj>dj,  {i.j)  e  Ek 


P{k,mo) 


where  1  is  the  set  of  jobs  to  be  processed  on  machine  k,  r,  is  the  release  time  (or  head)  of 
job  i,  equal  to  the  length  L{0,i)  of  a  longest  path  from  node  0  to  node  i  in  Z)(A/(i),  and 
is  the  delivery  time  (or  tail)  of  job  i,  equal  to  L(i,n)  —  d;,  where  L{iy7i)  is  the  length  of  a 
longest  path  from  node  i  to  node  n  in  D{Mo). 

On  the  other  hand,  the  one  machine  problem  with  DPC’s  for  the  same  machine  k  can  be 
stated  as 

minimize  t„ 

subject  to 

in  ~  ti  ^  dj  +  i  ^  I 

DPC{k,Mo)  U  >  r,  i  G  1 

ij-U  >  [i.j]  €  F 

V  {i,j)  G  Ek, 


where  F  is  the  set  of  precedence  arcs,  i.e.  pairs  with  delayed  precedence  constraints,  and 
L{i,j)  is  the  length  of  a  longest  path  from  i  to  j  in  the  D{Mo). 

It  is  a  well  known  fact  (see  [GJ79])  that  the  standard  one  machine  problem  P{k^Mo) 
is  NP-complete  in  the  strong  sense.  Since  the  problem  DPC{k,  Mq)  is  a  generalization  of 
P{k,Mo)  (it  becomes  P{k,Mo)  when  F  =  0),  it  is  also  NP-complete  in  the  strong  sense. 
The  question  that  we  want  to  address  now,  is  how  much  more  difficult  is  DPC{k,  Mo)  than 
P{k,Mo)l 

To  gain  some  insight  into  this  issue,  we  recall  two  well  known  facts  about  P{k,Mo)-  We 
call  the  Longest  Tail  rule  the  scheduling  of  jobs  in  order  of  decreasing  tail  length,  and  the 
Shortest  Head  rule  the  scheduling  of  jobs  in  order  of  increasing  head  length. 

Fact  1.  If  r,  =  Vj  for  all  i,  j  G  /,  then  the  Longest  Tail  rule  yields  an  optimal  schedule. 
If  qi  =  qj  for  all  i,  j  G  /,  then  the  Shortest  Head  rule  yields  an  optimal  schedule. 

Fact  2.  Fact  1  remains  true  if  preemption  (job  splitting)  is  allowed. 

Note  that  Longest  Tail  (Shortest  Head)  scheduling  requires  0(7?  log  n)  time. 

If  standard  (as  opposed  to  delayed)  precedence  constraints  are  imposed  on  the  one  ma¬ 
chine  problem,  i.e.  constraints  of  the  form  G  -  L  >  d,  for  some  (?,j)  G  F’t,  Facts  1  and  2 
remain  true,  provided  that  the  Longest  Tail  (Shortest  Head)  rule  is  applied  subject  to  the 


6 


precedence  constraints.  However,  in  the  presence  of  delayed  precedence  constraints,  both 
of  the  above  properties  break  down  and,  as  shown  by  the  next  two  theorems,  there  is  no 
polynomial-time  algorithm  even  for  the  more  restricted  case  where  r,  =  Vj  and  for 

all  i,  j  6  /,  unless  P  =  NP. 

For  the  next  two  Theorems  we  will  use  a  slightly  modified  version  of  the  DPC’s.  Namely, 
instead  of  requiring  that  the  start  time  of  job  i  precede  the  start  time  of  job  j  by  at  least 
L{i,j)  units,  we  will  require  that  the  completion  time  of  job  i  precede  the  start  time  of  job  j 
by  at  least  L'{i^j)  :=  L{i,j)  —  d,  units.  For  the  one  machine  scheduling  problem  with  DPC’s 
in  which  preemption  is  not  allowed  this  makes  no  difference,  since  the  completion  time  of 
job  i  is  Ci  :=  ti  +  d^,  hence  the  two  forms  of  the  DPC,  tj  —  U  >  L{i^j)  and  tj  —  C,  >  L'{i,j) 
are  equivalent.  But  in  the  case  when  preemption  is  allowed  this  is  no  longer  true,  since  a 
preempted  job  i  will  be  completed  later  than  its  starting  time  plus  its  processing  time;  so 
the  new  condition  yields  a  stronger  relaxation  of  the  one  machine  problem  with  DPC’s. 

Theorem  2.1  DPC{k,  Mq)  is  NP-hard  in  the  strong  sense  even  if  Vi  =  and  qi  =  q^  for 
all  i,j  €  I . 

Proof.  Our  proof  is  based  on  a  reduction  of  3-PARTITION.  The  problem  is  defined  as 
follows;  Given  3n  +  1  positive  integers  ai, . . . ,  a3„,  6,  with  bjA  <  a^  <  b/2  for  i  =  l,...,3n, 
and  a,  =  nb,  does  there  exist  a  partition  of  {!,...  ,3rr}  into  n  pairwise  disjoint  subsets 
5i , . . . ,  5n  such  that  ITigSh  Oi  =  6  for  /i  =  1 , . . . ,  n? 

We  will  show  that  for  every  instance  of  3-PARTITION  we  can  compute  an  instance  of 
DPC{k,  Mo)  and  a  value  2  in  polynomial  time,  such  that  the  instance  of  the  former  problem 
has  a  yes  answer  if  and  only  if  the  instance  of  the  latter  problem  has  a  schedule  of  length  at 
most  z.  The  Theorem  then  follows  from  the  known  fact  that  3-PARTITION  is  NP-complete 
in  the  strong  sense  [GJ79]. 

Given  an  instance  of  3-PARTITION  as  stated  above,  we  construct  the  following  instance 
of  DPC{k,  Mo).  There  are  4n  4-  1  jobs.  All  their  release  and  delivery  times  are  0.  The  first 
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3n  jobs  have  processing  times 


di  =  Cj  for  i  =  1, . . . ,  3n; 

there  are  no  precedence  constraints  between  these  jobs.  The  remaining  n  +  1  jobs  have 
processing  times 

di  =  b  for  2  =  3n  +  1 , . . . ,  4n  +  1 
and  also  are  subject  to  the  following: 

i  has  to  precede  i  +  1  with  L'{i,  f  +  1)  =  6  for  f  =  3ri  +  1, . . . ,  4n. 

Finally,  we  define  z  =  2nb  +  b.  Note  that  the  precedence  constraints  define  a  single  chain. 

Suppose  now  that  the  instance  of  3-PARTITION  is  a  yes  instance,  i.e.  there  exist  n 
pairwise  disjoint  subsets  Sh  with  for  h  =  1, . . . ,  n.  We  then  construct  a  schedule 

of  length  z  as  follows.  The  jobs  3n  +  1, . . .  ,4n  +  1  are  started  as  early  as  possible,  subject 
to  the  delayed  precedence  constraints,  i.e.  =  {2h  —  2)6  for  6  =  1, . . . ,  n  -f  1.  Note  that 
job  4n  -f  1  finishes  at  time  z.  This  leaves  n  idle  intervals  [(2A  —  1)6,266],  6  =  1, . . .  ,n,  of 
length  6  each,  in  which  we  process  the  jobs  i  with  i  €  5^  for  6  =  1, . . . ,  n. 

Conversely,  suppose  that  a  schedule  of  length  z  exists.  In  such  a  schedule  we  must  have 
tsh+h  =  (26  —  2)6  for  6  =  1, . . .  ,n  -(-  1:  none  of  these  jobs  can  start  earlier  in  view  of  the 
DPC’s,  and  none  of  them  can  start  later  in  view  of  the  schedule  length.  It  follows  that  all 
of  the  jobs  1, . . .  ,3n  are  processed  in  the  n  intervals  [(26  —  1)6,266],  6  =  1, . . .  ,n.  Now’ 
define  Sh  as  the  set  of  jobs  processed  in  the  6'^  of  these  intervals.  We  have  o-t  —  b  for 

6  =  1, . . .  ,n,  and  hence  we  have  a  yes  instance  of  3-PARTlTION.  □ 

Theorem  2.2  Theorem  2.1  remains  true  even  if  preemption  is  allowed. 

Proof.  We  again  give  a  reduction  of  3-PARTITION.  The  reader  should  verify  that,  w’hen 
preemption  is  allowed,  the  reduction  given  in  the  proof  of  Theorem  2.1  no  longer  w’orks:  a 
schedule  of  length  always  exists,  irrespective  of  the  answer  to  the  3-PARTITION  instance. 
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Given  an  instance  of  3-PARTITION,  we  construct  the  following  instance  of  the  preemp¬ 
tive  scheduling  problem  with  DPC’s,  with  release  and  delivery  times  all  equal  to  0.  There 
are  8n  1  jobs.  For  each  a,,  there  are  now  two  jobs,  i  and  3n  -t-  i,  with 


di  =  dsn+i  =  a,  for  i  =  1, . . . ,  3?; 


and  the  condition  that 

i  has  to  precede  3n  -f  i,  with  V{i,3ri  -f-  i)  =  2nb  —  6  for  z  =  1 . 37(. 

The  remaining  jobs  have  processing  times 

dt  =  b  for  z  =  6n  -f  1 , . . . ,  8n  -f  1 

and  are  subject  to  the  condition  that 

i  has  to  precede  z  +  1,  with  -f-  1)  =  6  for  z  =  6n  +  1, . . .  ,8zz. 

Finally,  we  define  r  =  4nb  +  b.  Note  that  the  precedence  constraints  define  a  collection  of 
chains  (one  chain  of  length  2n  —  1  and  3n  chains  of  length  1). 

Suppose  we  have  a  yes  instance  of  3-PARTITION.  A  schedule  of  length  .r  is  then  con¬ 
structed  as  follows.  We  start  the  jobs  6rz  -1-  1, . . . ,  8n  -|-  1  as  early  as  possible,  i.e.  ten+k  = 
{2h  —  2)b  for  h  =  1, . . . ,  2rz  -j-  1.  Again,  the  last  job  finishes  at  time  c.  We  are  left  with 
2rz  intervals  [(2/i  —  1)6, 2/z6]  of  length  6  each  for  the  jobs  l,...,6n.  For  k  =  l,....zz,  we 
process  the  jobs  z  with  i  G  Sh  in  [(2rz  -  1)6, 2n6],  and  their  successors  3zz  -|-  z  with  z  €  Sh  in 
[(2(tz  +  h)  —  1)6, 2(n  -f  6)6];  this  is  feasible  since  the  time  span  from  the  end  of  the  former 
interval  to  the  beginning  of  the  latter  is  equal  to  the  required  minimum  delay  of  2nb  -  b. 

Conversely,  suppose  that  we  have  a  schedule  of  length  r.  As  before,  we  must  have 
kn+h  =  (26  -  2)6  for  6  =  1, .  .  .  ,2n  -t-  1.  Hence  all  of  the  jobs  1, .  .  .  ,  (izz  are  processed  in 
the  intervals  [2(6  -  1)6.266],  6  =  l,...,2n.  Note  that  the  total  length  of  these  intervals 
equals  the  total  required  processing  time  of  these  jobs,  so  that  the  schedule  contains  no  idle 


time.  Also,  in  view  of  the  DPC’s,  the  jobs  must  be  processed  in  the  first  half  of 

the  schedule  and  their  successors  3n  +  1, . . .  ,6/j  in  the  second  half. 

Consider  the  first  interval,  [6.26],  and  let  be  the  set  of  jobs  started  and  completed  in 
this  interval.  Assume  that  a;  =  6  —  c  for  .some  c  >  0.  The  only  jobs  ‘.in  +  i  that  are 

available  for  processing  in  [(2?r  +  1)6,  (2n  +  2)6]  are  those  for  which  i  6  and  hence  that 
interval  contains  c  units  of  idle  time.  It  follows  that  our  assumption  is  incorrect  and  that 
=  b.  An  iterative  application  of  this  argument  leads  to  the  identification  of  sets 
5/i  with  IZjgSfe  di  =  6  for  A  =  2, . . . ,  n,  and  to  the  conclusion  that  we  have  a  yes  instance  of 
3-PARTITION.  □ 

As  mentioned  in  the  introduction,  the  standard  one  machine  problem,  although  NP- 
complete,  can  in  practice  be  solved  efficiently  for  fairly  large  sizes.  This  can  be  done  by  an 
algorithm  developed  by  Carlier  [C82],  based  on  an  analysis  of  the  properties  of  the  Longest 
Tail  Heuristic  (for  the  latter  see  also  Potts  [P80]).  The  Longest  Tail  Heuristic  for  the  standard 
one  machine  problem  schedules  the  jobs  sequentially,  choosing  at  each  step  the  job  with  the 
longest  tail  among  those  available  for  scheduling.  If  C  is  a  critical  path  in  the  digraph 
corresponding  to  the  resulting  schedule,  let  p  be  the  first  node  encountered  when  traversing 
C  backwards  starting  from  n  (not  counting  n),  let  c  be  the  first  node  (if  any)  encountered 
such  that  qc  <  qp,  and  let  J  be  the  set  of  nodes  traversed  (not  counting  n  and  c).  Potts' 
and  Carlier’s  analysis  established  that  (i)  if  there  is  no  node  c  with  the  above  property,  the 
schedule  at  hand  is  optimal;  and  (ii)  otherwise  in  any  optimal  schedule  job  c  comes  either 
before  or  after  all  the  jobs  in  J . 

Thus  Carlier’s  branch  and  bound  algorithm  applies  the  Longest  Tail  Heuristic  to  the 
problem  at  hand  at  any  given  node  of  the  search  tree.  If  the  resulting  schedule  is  optimal 
(i.e.  satisfies  (i)  above),  the  subproblem  corresponding  to  the  node  is  discarded;  othe'rwise 
it  is  replaced  by  two  new  subproblems,  obtained  from  the  parent  problem  by  im[)osing 
the  condition  that  job  c  has  to  precede  (first  descendant).  respectiv('ly  succcchI  (second 
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descendant)  all  jobs  in  J .  Carlier  also  uses  the  fact  that  for  any  I  C  N  the  value 

h{l)  :=  min  Vi  +  52*^'  +  ’ti”* 

‘e/  ^  .e/ 

is  a  valid  lower  bound  on  the  makespan,  and  applies  it  to  the  subproblems  generated  by 
branching,  with  /  =  J  for  the  first  and  7  =  J  U  {c}  for  the  second  descendant. 

The  efficiency  of  Carlier's  procedure  is  due  to  the  fact  that  in  the  absence  of  property 
(ii)  above,  all  the  situations  in  which  job  c  precedes  some  but  not  all  the  jobs  in  J  would 
have  to  be  considered. 

It  is  easy  to  see  that  in  the  presence  of  standard  (as  opposed  to  delayed)  precedence 
constraints,  properties  (i)  and  (ii)  still  hold.  Therefore  Carlier’s  branch  and  bound  algo¬ 
rithm  can  be  (and  has  been)  modified  to  accommodate  such  constraints  without  any  loss  in 
efficiency. 

This,  however,  is  not  the  case  when  delayed  precedence  constraints  are  imposed,  as  will 
be  seen  in  the  next  section. 


3  The  Longest  Tail  Heuristic  in  the  Presence  of  DPC’s, 
and  Its  Analysis 

We  will  now  state  a  m.odified  version  of  the  Longest  Tail  Procedure,  that  takes  into  account 
the  presence  of  DPC’s,  and  analyze  its  properties.  We  denote,  as  before,  by  1  the  set  of  jobs 
to  be  processed  on  the  given  machine.  Further,  we  introduce 

S  ;=  {i  ^  1  :  i  has  been  scheduled  } 

7r(i)  :=  {j  €  /  :  j  is  required  to  precede  i} 

(=  predecessors  of  i) 

<^{i)  '■=  {j  G  /  :  j  is  required  to  succeed  i} 

(=  succes.sors  of  i) 

The  sets  ;r(r)  and  a{i)  can  be  used  to  tighten  the  release  times  and  tails  by  capturing  the 
logical  implications  of  the  precedence  constraints.  Thus  the  numbers  can  be  updated 

as  follows; 


Tj  ;=  max{r,,  max{r,  -f  L(i,  j)}, 

<e7r(j) 


min  r,  +  ^  d,} 


(1) 


U 


and 


qj  max{qj,  max{I(j,i)  -  dj  +  d,  +  q,],  V  d,  +  min  q,}  (2) 

.eaU)  »erU) 

Longest  Tail  Heuristic  for  DPC’s 

Initialization.  Set  t  :=  0.,  S  :=  0,  r'  =  r,,  i  ^  1 . 

Iterative  Step.  Let  Q  be  the  set  of  unscheduled  jobs  whose  predecessors  have  all  been 
scheduled,  i.e. 

Q:={ieI\S:  7r(f)  C  S} 

Choose  k  such  that 

maxj^i  :  i  E  Q  and  r'-  <  t]  if  such  i  exists 

qk  ■=  < 

m&x{qi  :  i  E  Q  and  r'  :=  min{r'  :  i  €  Q})  otherwise. 

Set  tk  :=  maxjt,  r'f.},  S  :=  S  Li  {k} ,  t  tk  +  dk- 
For  all  j  G  cr{k),  set  r'  max{r',  4  +  L{k,j)]. 

If  /  \  5  ^  0,  repeat  the  Iterative  Step.  Otherwise,  define  the  set  of  delayed  jobs 

D  :=  {i  e  I  :  r'i>  n} 

and  stop. 

The  outcome  of  the  longest  tail  procedure  is  a  schedule  t  :=  (t,)  which  can  be  represented 
by  a  directed  graph  G{t)  with  one  or  several  longest  paths  (critical  paths)  from  0  to  n. 
Figure  3  shows  such  a  solution  digraph  for  a  problem  with  6  jobs.  The  critical  path  (0, 4,2,6, 7) 
is  drawn  in  heavy  lines.  The  component  of  <,  i.e.  the  start  time  of  the  job  (not  shown 
in  Figure  3),  is  equal  to  the  length  of  a  longest  path  from  0  to  node  i.  The  makespan  is  the 
length  of  a  longest  path  (critical  path)  from  0  to  n. 

For  the  standard  one  machine  problem,  it  is  a  well  known  fact  that  for  any  subset  K  C  /, 
the  number 

h{I\)  :=  minr,  +  V' d,  +  ming, 
leK  t€/\ 

leA 

is  a  lower  bound  on  the  makespan.  Furthermore,  the  strongest  such  bound,  max{//(A  )  : 
K  C  /},  is  known  to  be  equal  to  the  minimum  makespan  of  the  preemptive  version  of 
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Figure  3;  Digraph  of  a  schedule,  with  critical  path 

the  problem.  It  is  easy  to  see  that  these  bounds  are  also  valid  for  our  problem,  of  which 
the  standard  problem  is  a  relaxation.  In  particular,  if  C  :=  {(0,  ii),  (m,  ?2)i  •••)(«?' ^ 
critical  path  in  G{t)  and  Ic  •'=  then 

h{Ic)  :=  minr,  +  V]  d,  +  ming, 

is  a  lower  bound  on  the  makespan. 

Now  let  t  be  the  schedule  produced  by  the  Longest  Tail  Procedure,  G[t)  the  associated 
digraph,  and  C  ;=  {(0,  i\),  (ij,  12),  ■  ■  ■ ,  (ip,  n)}  a  critical  path  in  G(t),  with  Ic  ’■=  {?i, . . . ,  ip}- 
Further,  let  c  be  the  first  node  encountered  when  traversing  C  backwards  from  node  n  such 
that  qc  <  qip,  if  such  c  exists,  or  else  let  c  =  0,  with  tg  ;=  0.  Finally,  let  J  be  the  set  of  nodes 
of  C  between  n  and  c  (not  counting  n  and  c).  We  then  have  the  following 

Theorem  3.1  Suppose  >  max{<,,  i  €  J ,  and  the  segment  C(c,ip)  of  tli(  rnliral  path 
C  between  c  and  ip  contains  no  precedence  arc.  Then 
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(i)  If  c  =  0,  t  is  optimal; 


(ii)  If  c  >  Q,  then  in  any  schedule  bette7-  than  t.  job  c  either  precedes  or  succetds  all  jobs 
in  J. 

Proof,  (i)  Let  c  =  0.  By  assumption  r,  >  r,,  for  each  i  €  J.  hence  for  all  i  €  Ic  '-  hence 

r,,  =  min  r,.  Also,  since  the  arc  (0,  'A)  is  on  a  critical  path,  L,  =  •  Further.  >  q,j,  for 

ieic 

all  i  €  Ic  by  hypothesis,  hence  q,  =  min  q,.  Thus  the  lower  bound  h(Ic)  is  in  this  case 

t&lc 

+  E  d,  +  qtp.  But  this  is  precisely  the  length  of  C,  i.e.  the  makespan,  since  =  7-,, . 

»6/c 

+  d,^  ioT  k  =  I,. . .  ,p  -  1,  and  tn  =  Lp  +  d,^  +  q^^,.  Thus  /  is  optimal. 

(ii)  Suppose  now  that  c  >  0,  and  let  t’  be  a  schedule  in  which  c  comes  after  some  (but 
not  all)  jobs  in  J .  Since  the  total  processing  time  of  all  jobs  in  J  U  {c}  is  in  both  schedules 
the  same,  a  d  the  tail  of  the  last  job  in  f  is  at  least  as  long  as  qi^,  we  can  have  t'^  <  t^ 
only  if  <  tc  for  some  j,  €  J,  which  would  imply  r^,  <  tc-  But  >  tc  for  all  j  E  J  by 
hypothesis,  hence  t'^  >  tn,  i.e.  t'  is  not  oetter  than  t.  □ 

Notice  that  the  condition  >  max{Li,^c}  for  all  i  €  J,  always  satisfied  in  a  schedule 
produced  by  the  Longest  Tail  Heuristic  for  the  standard  one  machine  problem,  may  not  hold 
here  and  therefore  needs  to  be  explicitly  imposed.  Indeed,  although  the  updated  release 
times  r'  generated  during  the  procedure  always  satisfy  r'  >  max{L,  for  all  i  G  J  (which 
explains  why  job  i  was  not  chosen  before  jobs  i]  or  c),  one  may  still  have  r,  <  max{L, -fc}- 
Theorem  3.1  identifies  the  situation  in  which  the  same  efficient  branching  rule  used  in 
Carlier’s  algorithm  can  be  used  in  the  presence  of  DPC’s.  Next  we  address  the  situation  in 
which  the  assumptions  of  Theorem  3.1  do  not  hold. 

First  of  all,  note  that  if  a  critical  path  C  in  G(t)  contains  a  precedence  arc  (i,j)  such 
that  rj  =  tj,  then  {(0,y)}  U  C{j,ri)  is  a  critical  path  that  does  not  contain  {i,j)-  Similarly, 
if  C  contains  a  precedence  arc  {i,j)  such  that  q,  =  L{0,n)  —  ti  —  d,,  where  L(0,n)  is  the 
makespan  associated  with  /,  then  C’fO,  7)U{(7,  n)}  is  a  critical  path  not  containing  {i,j)-  Such 
precedence  arcs  are  of  no  consequence,  in  the  sense  that  a  critical  path  containing  such  an 
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arc  can  be  replaced  by  one  not  containing  it.  Also,  the  DPC's  discussed  here  are  not  the  only 
kind  of  precedence  constraints  that  can  occur.  Whenever  some  disjunction  can  be  arbitrated 
as  a  result  of  logical  tests  (see  the  next  section),  the  replacement  of  the  corresponding  edge 
with  a  directed  arc  creates  a  (nondelayed)  precedence  arc  with  a  constraint  of  the  type 
tj  —  ti  >  di.  These  nondelayed  precedence  arcs  are  also  of  no  consequence.  In  order  to 
distinguish  those  precedence  arcs  that  matter,  we  will  call  a  precedence  arc  {i,j)  tssential 
(with  respect  to  t)  if  L{i,j)  >  di,  <  tj  and  qi  <  Z/(0, n)  —  t,  —  d,. 

Proposition  3.2  If  every  critical  path  of  G{t)  contains  a  precedence  arc.  then  every  critical 
path  of  G{t)  contains  an  essential  precedence  arc. 

Proof.  Suppose  every  critical  path  of  G{t)  contains  a  precedence  arc,  but  there  is  a  critical 
path  C  that  contains  no  essential  precedence  arc.  Let  {i,j)  be  the  first  precedence  arc  of  C . 
Since  {i,j)  is  not  essential,  either  rj  =  tj  or  qi  =  1,(0,  n)  -  L  —  d,.  Let  C'  {(0,  j))  U  C(j.n) 
in  the  first  case,  and  C'  ;=  C(0,j)U  {(i,n)}  in  the  second.  In  either  case.  C'  is  a  critical  path 
containing  at  least  one  less  precedence  arcs  than  C.  Repeating  this  argument  for  as  long  as 
the  critical  path  at  hand  continues  to  have  a  precedence  arc  yields  a  critical  path  without 
any  precedence  arc,  a  contradiction.  □ 

From  now  on  we  focus  on  essential  precedence  arcs. 

The  reasoning  used  in  the  proof  of  Theorem  3.1  can  also  be  used  to  identify  other  situa¬ 
tions  in  which  certain  jobs  have  to  be  scheduled  before  (or  after)  other  jobs  in  any  schedule 
better  than  t.  Recognizing  such  situations  and  acting  on  them  can  in  turn  be  used  to  get 
rid  of  certain  essential  precedence  arcs. 

Proposition  3.3  Suppose  the  segment  C{j,n)  of  C.  where  j  =  q.  contains  no  precedence 
arc.  Let  K  :=  •  •  •  ,ip}  o.nd  suppose  qk  >  for  all  k  G  A  U  {j  },  and  rk  >  fj  for  all 

ke  I<\  cr{j). 

Then  in  any  schedule  better  than  t,  job  j  precedes  all  jobs  in  1\  . 
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Proof.  Let  t'  be  a  schedule  in  which  job  j  comes  after  some  or  all  of  the  jobs  in  A'.  Since  the 
total  processing  time  of  all  jobs  in  A  U  {j}  in  t'  cannot  be  less  than  what  it  is  in  t,  namelv 
^  d,,  and  the  tail  of  the  last  job  in  t'  is  by  hypothesis  at  least  etjual  to  .  we  can  have 

k^Ku{j} 

<  tn  only  if  <  tj  for  some  j.  €  A',  which  would  imply  <  tj.  But  according  to  our 
assumptions,  r*.  >  tj  for  all  k  G  h  \  u'(j),  and  for  k  6  crij)  k  cannot  precede  j  in  t'.  Hence 
t'n  ^  ^ni  i-e.  t'  is  not  better  than  t.  □ 

Now  suppose  the  segment  C(j,n)  with  the  properties  stated  in  Proposition  3.3  is  preceded 
on  the  critical  path  C  by  an  (essential)  precedence  arc  (Aj)  =  ?/).  Then  in  any  schedule 

better  than  t,  job  i  precedes  all  the  jobs  in  {j}  U  A'.  Hence  one  can  set 

cr{i)  (T{i)UK, 

7r{k)  :=  7r(A)U{z}.  k^K 

and  the  tail  of  job  i  can  be  increased  to 

qt  :=  -L(Ai)  +  dj  +  +  9ip  —  d, 

=  A(0,n)-A-d,,  ^ 

which  gives  rise  to  a  new  critical  path,  C(0,  i)U{(A«)},  that  does  not  contain  the  precedence 
arc  Thus  a  situation  like  this  can  be  used  to  render  an  essential  precedence  arc 

inessential. 

Next  we  consider  the  converse  of  the  above  situation. 

Proposition  3.4  Suppose  the  segment  C{0,i)  of  C,  where  i  =  A,  contains  no  precedence 
arc.  Let  K  :=  {A, .  .  . ,  A_i}  and  suppose  >  r,,  for  all  k  £  A'U  {i}  and  qk  >  1(0,  n)  —  t,  —  dj 
for  all  k  ^  K  \  7r(i). 

Then  in  any  schedule  better  than  t,  job  i  succeeds  all  jobs  in  K . 

Proof.  The  makespan  associated  with  t  is  <„  =  r,,  +  '^dh  +  A(^,n),  where  L{i.n)  — 

heK 

A(0,n)  —  A  is  the  length  of  the  critical  path  .segment  C{i.n}.  Now  let  t'  be  a  schedule  in 

which  job  i  precedes  some  or  all  of  the  jobs  in  A'.  Since  the  processing  time  of  all  jobs  in  A' 

in  t'  cannot  be  less  than  ^d/i  and  the  relea.se  time  of  each  job  in  A  U  {?  }  is  by  hypothesis 

heK 
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at  least  equal  to  r,, ,  it  follows  that  max  {t',^  d^}  >  ^  +  d,.  Therefore  we  can  have 

h&Ku{i} 

<  tn  only  if  the  part  of  the  schedule  t'  starting  after  completion  of  the  job  k  such  that 

t',  dk  :=  max  {t'f^  +  d^},  is  shorter  than  the  corresponding  part  of  the  schedule  t  (i.e.  the 
*r€/vu{i} 

part  of  t  starting  after  the  completion  of  job  i),  whose  length  is  L{i,n)  —  d,  —  L{0.n)  —  t,  —  d,. 

But  since  qt;  >  L{0,n)  —  —  d,  for  all  k  G  K  \  Tr{i)  by  assumption,  and  i  cannot  precede 

k  e  7r(i)  in  this  is  impossible  and  therefore  i.e.  is  not  better  than  ?/.  □ 

Similarly  to  the  previous  case,  if  the  segment  t'(0,?)  that  is  the  object  ol  Proposition  .3.4 

is  followed  on  the  critical  path  C  by  an  (essential)  precedence  arc  (ij)  =  (q.q+i).  then  in 

any  schedule  better  than  T  job  j  succeeds  all  the  jobs  in  A  U  {i}.  Hence  one  can  set 

:=  7r(j)UA'  .. 

a{k)  :=  <T(A;)U{j},  k^K 

and  the  release  time  of  job  j  can  be  increased  to 

r-j  :=  d^  +  L{iJ),  (6) 

h€K 

which  again  gives  rise  to  a  new  critical  path,  {(0,_;)}uC(j,  n),  not  containing  the  precedence 
arc  Thus  we  have  a  second  way  of  rendering  inessential  a  formerly  essential  precedence 

arc. 

4  Branch  and  Bound  for  the  One  Machine  Problem 
with  DPC’s 

We  implemented  a  branch-and-bound  procedure  that  w'orks  as  follow's. 

0.  Initialization.  The  original  problem  is  put  on  the  list  of  active  subproblems,  after 
updating  the  release  times  and  tails  according  to  expressions  (1 ).  (2)  of  section  .3,  initializing 
the  upper  bound  at  f  :=  oc  and  calculating  a  lower  bound  h  :=  max  h.(Q),  where 

h(Q)  :=  minr,  +  Y"  d,  -f-  minq,  (7) 

.eg  ^  .eg 

It  is  known  (see  [C82])  that  h(Q)  can  be  calculated  in  O(nlogn)  time  by  solving  the 
preemptive  version  of  the  standard  one  machine  problem. 
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1.  Subproblem  Selection.  The  problem  with  the  weakest  (i.e.  smallest)  lower  bound  is 
selected  for  processing  next,  and  is  removed  from  the  list. 

2.  Longest  Tail  Heuristic.  The  selected  problem  is  processed  by  the  longest  tail  heuristic 
of  section  3  to  obtain  a  schedule  t.  If  <  /,  the  upper  bound  is  updated  to  /:=/„. 

3.  Postprocessing.  Propositions  3.3  and  3.4  are  applied:  the  predecessor  and  successor 
sets,  and  the  release  times  and  tails,  are  updated  by  using  expressions  3,  4.  5  and  6.  If  any 
tail  has  been  changed,  go  to  2;  otherwise  go  to  4. 

4-  Branching.  We  search  for  a  critical  path  satisfying  the  conditions  of  Theorem  3.1,  i.e. 
one  that  contains  no  precedence  arc  in  the  .segment  C(c,n)  and  for  which  r,  >  max{t,, , 
for  ail  i  ^  J  D  D.  If  found,  then  either  c  —  0  and  the  schedule  t  is  optimal  for  the  given 
subproblem,  hence  the  latter  is  discarded  and  we  go  back  to  step  1;  or  else  c  >  0  and  we 
apply  the  strong  branching  rule  to  be  described  below.  If  the  search  fails,  we  repeat  it  on  the 
reverse  problem.  This  is  obtained  by  reversing  the  direction  of  every  arc,  interchanging  the 
release  times  with  the  tails,  and  reversing  the  precedence  constraints,  with  the  understanding 
that  if  the  delay  between  i  and  j  in  the  original  problem  is  L{i,j),  then  the  delay  between  j 
and  i  in  the  reverse  problem  is  L{j,i)  =  L{i,j)  —  di  +  dj.  The  reverse  problem  is  easily  seen 
to  be  equivalent  to  the  original  one;  thus  if  the  search  is  successful  on  the  reverse  problem, 
we  proceed  in  the  same  way  as  in  case  of  success  on  the  original  problem. 

If  the  search  fails  on  both  problems,  i.e.  either  for  every  critical  path  C  that  contains  no 
precedence  arc  in  the  segment  C{c,n)  we  have  r,  <  max{<,,  ,tc)  for  some  i  ^  J  H  D  (case  1), 
or  else  every  critical  path  C  contains  some  essential  precedence  arc  in  the  segment  C{c,ri) 
(case  2),  we  apply  the  weak  branching  rule  (see  below). 

a.  Strong  branching  rule.  We  replace  the  current  subproblem  by  two  new  ones;  a  left 
subproblem,  in  which  c  is  required  to  precede  all  jobs  in  J,  i.e. 

<t(c)  :=  o-(c)U./,  . 

7r{k)  :=  7r(yb)U{c},  k  £  J: 
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(9) 


and  a  right  subproblem,  in  which  c  is  required  to  succeed  all  jobs  in  7,  i.e. 

7r(c)  :=  7r(c)  U  J 

a{k)  ;=  <T{k)U{c},  k  G  J. 

The  strong  branching  rule  is  preceded  by  the  logical  tests  proposed  by  Carlier  [C82]  for 

the  standard  one  machine  problem,  which  obviously  preserve  their  validity  in  the  presence 

of  DPC’s  for  the  situation  when  the  strong  branching  rule  applies. 

Let  K  :=  {k  :  k  Ei  I  \  J  U  {c}  and  dk  >  f  —  h{J)}. 


Test  1.  \i  h  ^  K  and 

^  di  +  >  / 

(10) 

ieJ 

then  job  k  has  to  succeed  all 

the  jobs  in  J;  so  we  set 

7r(A:)  :=  ^(A:)  U  J 

a{i)  <j(f)U{Ar),  z  E  7. 

(11) 

Test  2.  If  A:  G  K  and 

min  r,  +  V  7i  +  +  gfc  >  / 

(12) 

then  job  k  has  to  precede  all 

the  jobs  in  7;  so  we  set 

cr(A:)  :=  a{k)Dj 

7r(z)  :=  7r(z)U{A;},  z  €  7. 

(13) 

b.  Weak  branching  rule.  Select  a  pair  of  jobs  {i,j)  as  follows.  In  case  1,  let  C  be 
one  of  the  critical  paths  without  precedence  arcs  in  C{c,n).  Let  j  =  c  c  ^  0,  j  =  i\ 
otherwise,  and  let  i  be  the  first  node  encountered  when  traversing  the  segment  C(c,  n), 
such  that  ri  <  max{i„,ic}-  In  case  2,  let  (7  be  a  critical  path  containing  a  minimum 
number  of  precedence  arcs  in  C(c,  n),  and  let  {k,l)  be  the  one  nearest  to  n.  Then,  {k,l)  is 
essential  (since  it  was  not  eliminated  by  postprocessing),  and  rm  <  ti  for  at  least  one  m  ^  J 
(otherwise  Proposition  3.3  would  apply).  Let  j  =  /  and  let  i  be  the  first  node  encountered 
when  traversing  the  segment  C(l,n)  such  that  r,  <  ti. 

In  both  cases,  we  replace  the  current  subproblem  by  two  new  ones:  a  left  subproblem,  in 
which  i  is  required  to  precede  j,  i.e. 

aii)  :=  o-(t)U{i}  .... 

7r(i)  :=  7r(j)U{?;}, 
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and  a  right  subproblem,  in  which  i  is  required  to  succeed  j,  i.e. 

7r(2)  :=  7r(z)U{j} 

o-(i)  := 

5.  Updating  and  Lower  Bounding 

For  each  of  the  two  new  subproblems  created,  we  update  the  release  times  and  tails  by 
using  expressions  (1)  and  (2).  In  fact,  we  use  a  strengthened  version  of  (1)  and  (2),  defined 
as  follows: 

Let  |x(j)|  =  m.  Let  Ho  :=  7r(j), 

rh{\)  ■-  min{r/^  :  h  e  Hq},  Hi  :=  Ho  \ 

and  for  A:  =  2, ... ,  m,  define  recursively 

rh{k)  ■=  min{rb  :  h  e  /ffc-i},  H^  :=  \  {h{k)}. 

Then  (1)  can  be  strengthened  to 

(F)  r,  :=  maxfr,,  maxfri  +  ?')},  max  {minr,  +  y^  di]]. 

^  i€H3)  '  fc€{o,i,...,m}Lew* 

Similarly,  if  lcr(j)|  =  s,  let  A'o  =  cr{j), 

:=  m\n{qf,  :  h  e  Kq},  A'l  :=  A'o  \  {/i(l)}, 

and  for  /  =  2, . . . ,  s,  define 

qh{i)  :=  minj^,.  :  h  e  A';_i},  A'/  :=  A';_i  \  {/i(/)}. 

Then  (2)  can  be  strengthened  to 

(2')  qj  :=ma.x{qj,  max{<7,  +  A(i,  j)  +  d,  -  d,},  max  {  V  d,  +  min^,}}, 

leKi 

Next  we  calculate  the  lower  bound  h{Q)  of  expression  (7)  by  solving  the  preemptive 
version  of  the  standard  one  machine  problem.  Finally,  we  put  each  of  the  two  subproblems 
on  the  list  of  active  problems  and  go  to  1. 
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5  Computational  Results 


The  branch  and  bound  method  of  section  4  was  implemented  in  C  and  tested  on  a  number 
of  randomly  generated  one  machine  problems  with  DPC’s.  It  was  then  embedded  into  the 
Shifting  Bottleneck  Procedure  in  place  ot  the  subroutine  for  solving  one  machine  problems. 
The  resulting  Modified  Shifting  Bottleneck  Procedure  (SB3)  was  tested  on  a  number  of  job 
shop  scheduling  problems  from  the  literature,  and  compared  with  other  procedures. 

We  first  discuss  the  computational  testing  of  the  algorithm  of  section  4  by  itself.  The  test 
problems  were  randomly  generated  in  the  same  way  as  those  of  [C82],  except  for  the  DPC's 
that  were  added.  We  set  ^  {50, 100},  r^ax  =  9inax  =  where  n  €  {20,50, 100} 

and  k  €  {10,15,20}.  Here  n  is  the  number  of  jobs  and  k  is  Carlier’s  “interval  multiplier.” 
The  numbers  ri,  qi,  d,  were  randomly  drawn  from  uniform  distributions  between  1  and  rmax, 
9tnax5  and  d^ax,  respectively.  A  precedence  constraint  was  generated  between  jobs  i  and  j 
with  a  probability  of  pij  €  {0,0.02,0.04, . . .  ,0.20},  and  the  delays  L(i,j)  were  generated  as 
follows:  for  each  {i,j)  €  F,  a  number  f(f,y)  was  drawn  from  a  uniform  distribution  over  the 
interval  [1,  ^n/:dmax];  and  L{i,j)  was  set  to  l{ij)  if  >  di,  and  to  +  di  otherwise. 
For  every  combination  of  and  fc,  40  instances  were  generated.  The  results  are  shown 
in  Tables  1-6.  As  one  can  see  from  these  tables,  the  number  of  search  tree  nodes  rarely 
exceeds  the  number  of  jobs,  which  is  in  contrast  with  the  typical  behavior  of  branch  and 
bound  algorithms  on  hard  problems,  but  is  comparable  to  (though  not  quite  as  good  as)  the 
performance  of  Carlier’s  algorithm  on  the  standard  one  machine  problem.  The  remarkable 
fact  is  that  typically  strong  branching  occurs  more  frequently  than  weak  branchings  and 
every  strong  branching  discards  a  large  number  of  potential  nodes  of  the  search  tree. 

The  main  test  of  the  usefulness  of  the  algorithm  for  the  one  machine  problem  with 
delayed  precedence  constraints  is  the  effect  on  the  performance  of  the  Shifting  Bottleneck 
Procedure  of  replacing  the  standard  one  machine  problem  by  the  one  machine  problem  w'ith 
DPC’s  as  the  tool  for  identifying  bottleneck  machines  and  scheduling  them.  This  modified 
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Table  1.  n  =  20,  <imax  =  50 


Density 

of 

DPC’s 

(%) 


Search 


Tree  SPARC330 

Nodes  Seconds 

Avg.  (Max.)  Avg.  (Max.) 


Proportion 
of  Strong 
Branchings 
(%) 


2.13  (  45) 
5.98  (  72) 
6.08  (  33) 

39.50  (  485) 

29.50  (  198) 

80.85  (1,818) 
16.38  (  100) 
29.75  (  201) 
16.58  (  134) 
20.20  (  107) 
16.55  (  126) 


5.53  (  21) 
8.88  (  25) 
14.35  (  135) 
14.60  (  124) 

7.85  (  18) 

8.38  (  24) 
10.75  (  66) 
6.30  (  18) 

6.63  (  23) 
5.13  (  22) 
4.90  (  18) 


11.93  (  168) 
6.45  (  15) 

6.33  (  22) 
5.28  (  16) 

3.63  (  8) 

4.00  (  13) 

4.18  (  14) 

3.85  (  24) 
3.30  (  35) 
3.60  (  11) 

3.00  (  18) 


0.008  (  0.183) 
0.027  (  0.400) 
0.028  (  0.217) 
0.293  (  4.500) 
0.220  (  1.567) 
0.587  (14.799) 
0.102  (  0.700) 
0.212  (  1.817) 
0.115  (  1.317) 
0.147  (  0.900) 
0.111  (  1.200) 


0.023  (  0.083) 
0.035  (  0.183) 
0.071  (  0.867) 
0.088  (  1.083) 
0.037  (  0.100) 
0.047  (  0.200) 
0.057  (  0.400) 
0.031  (  0.117) 
0.038  (  0.150) 
0.031  (  0.133) 
0.300  (  0.117) 


0.061  (  1.183) 
0.025  (  0.100) 
0.024  (  0.133) 
0.024  (  0.133) 
0.016  (  0.033) 
0.016  (  0.050) 
0.022  (  0.117) 
0.024  (  0.150) 
0.019  (  0.200) 
0.021  (  0.067) 
0.019  (  0.100) 


100.0 

81.7 


58.6 

50.4 
57.3 
51.1 

41.5 
50.8 

51.7 


100.0 

86.9 

65.5 

60.6 
86.4 

73.8 
67.1 

84.3 

75.9 

67.3 

65.4 


94.6 

100.0 

83.8 

89.1 

79.8 
90.3 

89.2 

86.2 


Table  2.  n  =  20,  =  100 


k 

Density 

of 

DPC’s 

(%) 

Search 

Tree 

Nodes 

Avg.  (Max.) 

CPU 

SPARC330 
Seconds 
Avg.  (Max.) 

1.43  (  9) 

0.006  (  0.0.50) 

15.28  (  239) 

0.091  (  1.250) 

8.15  (  73) 

0.043  (  0.600) 

6 

70.50  (  981) 

0.510  (  7.683) 

10 

8 

34.85  (  484) 

0.276  (  4.666) 

10 

46.55  (  459) 

0.342  (  4.700) 

12 

43.10  (  424) 

0.332  (  4.483) 

14 

32.65  (  324) 

0.215  (  2.550) 

16 

52.80  (1273) 

0.465  (12.450) 

18 

13.20  (  74) 

0.085  (  0.750) 

20 

11.80  (  126) 

0.080  (  1.183) 

0 

5.65  (  24) 

0.020  (  0.100) 

2 

15.18  (  157) 

0.081  (  1.183) 

4 

15.55  (  185) 

0.078  (  1.217) 

6 

13.68  (  47) 

0.074  (  0.350) 

15 

8 

9.28  (  49) 

0.048  (  0.417) 

10 

7.83  (  32) 

0.040  (  0.283) 

12 

9.93  (  64) 

0.050  (  0.400) 

14 

7.75  (  31) 

0.043  (  0.267) 

16 

6.08  (  21) 

0.034  (  0.150) 

18 

5.10  (  25) 

0.031  (  0.183) 

20 

5.30  (  18) 

0.031  (  0.100) 

0 

11.40  (  154) 

0.055  (  1.067) 

2 

6.43  (  17) 

0.023  (  0.067) 

4 

5.63  (  15) 

0.020  (  0.050) 

6 

5.30  (  18) 

0.024  (  0.117) 

20 

8 

3.78  (  9) 

0.017  (  0.033) 

10 

4.05  (  13) 

0.016  (  0.050) 

12 

4.20  (  15) 

0.021  (  0.117) 

14 

3.95  (  19) 

0.024  (  0.150) 

16 

2.63  (  15) 

0.017  (  0.100) 

18 

3.63  (  11) 

0.019  (  0.067) 

20 

2.90  (  15) 

0.020  (  0.067) 

23 


Table  3.  n  =  50,  c?„iax  =  50 


Density 

of 

DPC’s 

(%) 

Search 

Tree 

Nodes 

Avg.  (Max.) 

CPU 

SPARC330 
Seconds 
Avg.  (Max.) 

Proportion 
of  Strong 
Branchings 
(%) 

0 

1.03  ( 

2) 

0.014  (  0.050) 

100.0 

2 

34.60  ( 

163) 

0.587  (  3.117) 

54.8 

4 

36.48  ( 

259) 

0.646  (  5.550) 

47.9 

6 

18.40  ( 

58) 

0.333  (  1.067) 

56.5 

10 

8 

13.95  ( 

85) 

0.291  (  1.950) 

54.2 

10 

83.30  (1,752) 

2.449  (50.398) 

45.1 

12 

15.43  ( 

167) 

0.411  (  4.266) 

43.4 

14 

9.18  ( 

49) 

0.265  (  1.217) 

55.3 

16 

6.85  ( 

16) 

0.222  (  0.450) 

62.5 

18 

8.28  ( 

112) 

0.264  (  2.983) 

46.0 

20 

5.03  ( 

13) 

0.193  (  0.433) 

67.0 

32.90  ( 

798) 

0.980  (31.132) 

100.0 

2 

22.83  ( 

141) 

0.365  (  3.300) 

75.4 

4 

11.55  ( 

44) 

0.188  (  0.767) 

77.6 

6 

8.15  ( 

42) 

0.150  (  0.800) 

66.9 

15 

8 

7.10  ( 

33) 

0.156  (  0.667) 

70.5 

10 

4.83  ( 

12) 

0.125  (  0.317) 

78.6 

12 

5.25  ( 

24) 

0.152  (  0.817) 

69.4 

14 

3.60  ( 

12) 

0.123  (  0.317) 

87.0 

16 

4.88  ( 

18) 

0.172  (  0.500) 

64.4 

18 

4.18  ( 

23) 

0.168  (  1.017) 

75.8 

20 

4.05  ( 

16) 

0.156  (  0.450) 

75.0 

0 

18.58  ( 

39) 

0.238  (  0.533) 

100.0 

2 

6.98  ( 

19) 

0.103  (  0.300) 

90.4 

4 

5.28  { 

26) 

0.090  (  0.550) 

81.9 

6 

3.53  ( 

13) 

0.072  (  0.300) 

90.9 

20 

8 

4.05  ( 

9) 

0.095  (  0.200) 

81.5 

10 

3.43  ( 

10) 

0.096  (  0.217) 

84.8 

12 

2.98  ( 

10) 

0.100  (  0.250) 

93.2 

14 

3.10  ( 

11) 

0.112  (  0.317) 

86.4 

16 

2.83  ( 

8) 

0.112  (  0.233) 

81.1 

18 

2.45  ( 

10) 

0.116  (  0.267) 

82.0 

20 

2.50  ( 

7) 

0.119  (  0.267) 

87.2 

24 


Table  4.  n  =  50,  t/max  =  100 


1 

Density 

of 

DPC’s 

(%) 

Search 

Tree 

Nodes 

Avg.  (Max.) 

CPU 

SPARC330 

Seconds 

Avg.  (Max.) 

Proportion 
of  Strong 
Branchings 
(%) 

0 

1.05  (  2) 

0.019  (  0.067) 

100.0 

2 

34.03  (  123) 

0.578  (  2.383) 

53.6 

4 

39.15  (  179) 

0.694  (  3.567) 

45.2 

6 

99.53  (1866) 

2.814  (48.465) 

45.3 

10 

8 

36.23  (  899) 

0.860  (23.032) 

43.7 

10 

11.95  (  47) 

0.267  (  1.033) 

67.2 

12 

8.95  (  31) 

0.241  (  0.800) 

66.3 

14 

9.90  (  32) 

0.267  (  0.833) 

53.5 

16 

12.33  (  137) 

0.440  (  6.433) 

53.3 

18 

7.78  (  31) 

0.251  (  0.933) 

57.4 

20 

5.60  (  28) 

0.216  (  0.817) 

71.3 

0 

24.60  (  45) 

0.338  (  0.617) 

100.0 

2 

17.60  (  58) 

0.260  (  1.017) 

73.4 

4 

8.20  (  23) 

0.131  (  0.417) 

84.4 

6 

6.63  (  24) 

0.134  (  0.483) 

72.5 

15 

8 

5.95  (  27) 

0.134  (  0.500) 

69.2 

10 

4.48  (  18) 

0.118  (  0.417) 

79.3 

12 

3.60  (  11) 

0.112  (  0.483) 

84.1 

14 

4.38  (  12) 

0.155  (  0.667) 

85.1 

16 

3.75  (  17) 

0.134  (  0.417) 

82.4 

18 

3.43  (  15) 

0.143  (  0.467) 

82.1 

20 

3.83  (  11) 

0.163  (  0.367) 

74.7 

0 

19.18  (  35) 

0.240  (  0.450) 

100.0 

2 

9.18  (  19) 

0.123  (  0.267) 

90.8 

4 

4.65  (  14) 

0.074  (  0.233) 

84.3 

6 

3.55  (  11) 

0.071  (  0.167) 

82.8 

20 

8 

3.98  (  10) 

0.094  (  0.250) 

92.1 

10 

3.18  (  11) 

0.093  (  0.217) 

86.7 

12 

3.25  (  10) 

0.101  (  0.267) 

83.3 

14 

2.75  (  8) 

0.107  (  0.233) 

85.5 

16 

2.95  (  12) 

0.122  (  0.450) 

78.7 

18 

2.58  (  8) 

0.120  (  0.283) 

88.5 

20 

2.45  (  7) 

0.116  (  0.233) 

84.4 
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Table  5.  n  =  100,  c?max  =  50 


Density 

of 

DPC’s 

(%) 

Search 

Tree 

Nodes 

Avg.  (Max.) 

CPU 

SPARC.330 
Seconds 
Avg.  (Max.) 

Proportion 
of  Strong 
Branchings 
(%) 

0 

1.00  (  1) 

0.047  (0,050) 

100.0 

2 

48.65  (122) 

2.755  (7.566) 

44.4 

4 

15.00  (  58) 

0.955  (3.583) 

51.7 

6 

11.08  (  25) 

0.876  (1.800) 

54.6 

10 

8 

8.30  (  23) 

0.883  (2.200) 

62.0 

10 

7.90  (  19) 

0.955  (2.083) 

59.1 

12 

6.98  (  29) 

0.931  (3.667) 

63.2 

14 

6.15  (  18) 

0.907  (2.150) 

61.6 

16 

5.30  (  19) 

0.840  (2.217) 

64.7 

18 

4.23  (  14) 

0.795  (2.317) 

72.4 

20 

5.00  (  17) 

0.899  (2.417) 

63.4 

0 

42.50  (  93) 

1.920  (4.200) 

100.0 

2 

12.58  (  33) 

0.608  (1.650) 

76.8 

4 

6.78  (  35) 

0.450  (2.250) 

71.5 

6 

4.58  (  14) 

0.431  (0.983) 

78.9 

15 

8 

4.65  (  15) 

0.548  (1.667) 

92.1 

10 

4.43  (  16) 

0.615  (1.583) 

70.0 

12 

4.00  (  12) 

0.638  (1.283) 

83.5 

14 

3.76  (  11) 

0.655  (1.300) 

73.7 

16 

3.43  (  9) 

0.654  (1.200) 

89.2 

18 

3.13  (  8) 

0.671  (1.317) 
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20 

2.85  (  12) 

0.658  (1.650) 
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0 

31.70  (  55) 

1.314  (2.333) 

100.0 

2 

5.00  (  14) 

0.248  (0.717) 

92.6 

4 

4.43  (  19) 

0.324  (1.083) 

73.6 

6 

4.18  (  18) 

0.398  (1.317) 

76.5 

20 

8 
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2.35  (  5) 
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87.5 

18 

2.35  (  7) 

0.560  (0.950) 

88.4 

20 

2.08  (  6) 

0.570  (1.017) 

75.6 
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Table  6.  n  =  100,  t^max  =  100 


Density 

of 

DPC’s 

(%) 

Search 

Tree 

Nodes 

Avg.  (Max.) 

CPU 

SPARC330 
Seconds 
Avg.  (Max.) 

Proportion 
of  Strong 
Branchings 
(%) 

0 

1.00  (  1) 

0.049  (  0.133) 
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2 

49.95  (129) 

2.721  (  7.950) 

54.9 

4 
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60.8 

6 
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0.787  (  2.283) 

71.3 

10 
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8.43  (  29) 

0.908  (  3.467) 

54.3 

10 
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1.034  (  2.433) 

56.9 

12 

6.63  (  17) 

0.905  (  1.967) 

63.9 

14 

7.25  (  19) 

1.026  (  2.533) 

60.8 

16 

5.58  (  31) 

0.985  (  6.283) 

61.5 

18 

6.20  (  12) 

1.017  (  1.833) 

57.4 

20 

5.93  (  59) 

1.109  (10.916) 

63.8 

0 

55.20  (169) 

2.593  (  9.816) 

2 

10.18  (  29) 

0.481  (  1.350) 

92.3 

4 

7.85  (  27) 

0.539  (  1.667) 

76.9 

6 

6.43  (  20) 

0.595  (  2.083) 

68.7 

15 

8 

4.75  (  21) 

0.580  (  2.117) 

70.1 

10 

4.43  (  16) 

0.614  (  1.667) 

74.2 

12 

4.03  (  15) 

0.662  (  1.983) 

71.4 

14 

3.68  (  15) 

0.654  (  1.867) 

69.4 

16 

3.25  (  15) 

0.665  (  2.117) 

78.5 

18 

2.15  (  7) 

0.552  (  1.183) 

86.5 

20 

3.20  (  7) 

0.688  (  1.167) 

87.3 

0 

31.73  (  71) 

1.313  (  3.217) 

100.0 

2 

5.48  (  19) 

0.274  (  1.033) 

83.7 

4 

4.60  (  18) 

0.325  (  0.917) 

75.9 

6 

3.90  (  18) 

0.381  (  1.333) 

76.6 

20 

8 

2.83  (  7) 

0.395  (  0.883) 

90.4 

10 

2.95  (  10) 

0.483  (  1.100) 

77.2 

12 

4.08  (  17) 

0.663  (  2.083) 

73.2 

14 

2.65  (  9) 

0.523  (  1.050) 

87.0 

16 
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0.589  (  1.233) 

83.0 
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93.3 

20 

2.08  (  6) 

0.562  (  1.033) 

79.5 
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SB  Procedure,  which  we  cal!  SB3,  differs  from  SBl,  the  original  SB  Procedure  of  [ABZ88]. 
in  several  minor  details  besides  the  major  differences  of  using  the  one  machine  problem  with 
DPC’s,  and  the  number  of  reoptimization  cycles.  We  limit  the  number  of  reoptimization 
cycles  to  at  most  six. 

We  also  implemented  a  slightly  different  version  of  the  new  procedure  which  we  call 
SB4.  In  SB3,  upon  completion  of  the  local  reoptimization  procedure  for  a  given  set  A/o  of 
machines  already  sequenced,  the  procedure  is  repeated  after  temporarily  removing  from  the 
problem  the  last  a  non-critical  machines  (see  [ABZ88]  for  details).  In  SB4  after  applying 
this  local  reoptimization  as  described  abov^e,  we  apply  the  procedure  a  second  time  in  the 
reverse  order,  i.e.  by  first  temporarily  removing  from  the  problem  the  last  a  non-critical 
machines,  and  then  applying  again  the  reoptimization  step  with  all  the  machines  present. 
As  it  can  be  expected,  SB4  is  computationally  more  expensive  than  SB3  (roughly  by  a  factor 
of  2),  but  it  often  finds  a  better  solution. 

The  Procedures  SB3  and  SB4  were  tested  on  the  40  problems  generated  by  Lawrence 
[L84]  and  compared  to  SBl,  as  well  as  to  the  modified  SB  procedure  of  Dauziere- Peres 
and  Lasserre  [DL90],  which  also  uses  the  one  machine  problem  with  DPC's.  but  solves  it 
heuristically  rather  than  to  optimality.  We  also  compare  our  procedures  to  Applegate  and 
Cook’s  [AC91]  implementation  of  the  shifting  bottleneck  procedure.  The  results  are  shown 
in  Tables  7-9,  where  m  and  n  denote  the  number  of  machines  and  jobs,  respectively,  DL 
stands  for  the  procedure  of  Dauziere- Lasserre,  and  AC  for  that  of  Applegate  and  Cook.  The 
data  for  the  DL  procedure  were  kindly  communicated  to  us  by  Dauziere- Peres  [D91].  The 
data  for  the  AC  column  were  obtained  by  us  with  the  code  of  Applegate  and  Cook  [.'\C9l], 
kindly  provided  to  us  by  David  Applegate.  The  optimal  makespans  are  from  [AC91]. 

In  addition,  we  .solved  the  problems  FT1-FT3  generated  by  Fisher  and  Thompson  [FT63]. 
ABZ5-ABZ6  generated  by  Adams,  Balas  and  Zawack  for  [ABZ88].  and  the  problem.'^  ORBl- 
ORB5  generated  by  Applegate  and  Cook  for  [AC91].  We  ran  SBl,  SB3.  SB4  and  .4C  on 
these  problems,  with  the  outcome  shown  in  Table  10  and  Table  11.  Table  12  shows  the 
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Excluding  problems  ABZb-ORB 


number  of  times  each  procedure  ranked  first,  counting  and  not  counting  ties.  Since  some 
problem  sets  are  easy  for  all  three  procedures,  in  Table  13  we  show  the  same  data  for  the 
problems  2-4,  16-30,  36-40,  ABZ5-ORB5,  FT1-FT3,  that  remain  after  eliminating  the  easy 
sets. 

We  conclude  that  the  Modified  Shifting  Bottleneck  Procedure  (SB3).  and  its  variant 
(SB4)  find  consistently  better  schedules  than  SBl,  DL  or  AC,  at  a  computational  cost 
somewhat  higher  than  that  of  SBl  and  AC,  but  not  higher  than  that  of  DL. 
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